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INTRODUCTION 



Fatigue fractures occur at points of stress concentration, 
which are normally regions of complex stress distribution with 
large gradients. The stress field is geometry dependent and 
will only be reproduced in a laboratory sample if geometric 
similitude is preserved. Matching sample to part must also 
include similitude in the two major phases of fatigue; initia- 
tion and propagation and not merely a m.atching of stress con- 
centration factors. In fact if the fatigue process were well 
understood, it would seem that a knowledge of the local 
conditions, coupled with material property information, would 
be all that is required to predict fatigue life. 

This has not been possible so far, but this study has been 
undertaken to see how much success can be achieved toward this 
end via a detailed study of conditions surrounding the notch 
tip. To have any chance of success whatsoever, the local stress 
and strain conditions at the notch must be determined accu- 
rately. One of the means available for accomplishing this 
is experiment, and this approach was used successfully. 



An evaluation of how well one notch specimen models 
another can be studied via similitude theory. The strains in 
the vicinity of a notch depend upon the following variables : 



SIMILITUDE 




( 1 ) 
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= strain components 
= Applied load 
= Sheet width 
= Sheet thickness 
= Radius of curvature of notch 
= Depth of notch 
= Tensile Modulus 
= Poisson's Ratio 
By dimensional analysis, the functional relationship can be 
written in terms of groupings of variables that are dimen- 
sionless . 
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From similitude theory, if in two different geometries 
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GEO\!ETRIC SIMILARITY 



If the material of the part and the model are the same, 
the same strains are produced only if they are geometrically 
similar. It is not sufficient to match only the stress at 
the notch tip, but complete similarity must be maintained. 

If the complete strain fields are not the same, the plastic 
zones will be different, and the neighborhood where fatigue 
action is taking place has different conditions. Thus, 
similarity m.ust be preserved when employing experimental 
fatigue results from one specimen to predict fatigue failure 
in another. This is not generally done in current practice. 
Since it is not convenient, nor even practical, to run tests 
on a geometrically simiilar specimen each tim.e a part is an- 
alyzed for fatigue, a different approach is attempted to 
develop a more general theory where knowledge of local 
conditions at the notch tip is sufficient, along with 
material properties, to predict fatigue in any specimen. 

SIMILITUDE IN UNIAXIAL TENSION 

Because much of the fatigue data available has been 
collected on uniaxial tension specimens and because some- 
thing resembling uniaxial tension prevails at notch tips, 
it is reasonable to pursue this similarity, which is shown 
in Figure 1. 
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Figure 1. Uniaxial Model of Stress Condition at 
a Notch Tip. 

For an elastic material 

- ■(•I , ^y> S (3) 

If the same material is used for the model and the specimen, 
the strain field will be modeled to the degree that the 
stress field is modeled. For instance, in the region of 
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the notch tip in plane stress. 
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% 



L 

£ 




( 4 ) 



To keep the departure from a uniaxial stress state less than 
3%, which would be in the range of experimental measurement 
error, the transverse stress must be small. 

(T^ < 0,05 (5) 

For a typical value of Poisson's Ratio for aluminum, = 0.33 
and 

6"^ i O. O \ CTy (6) 

In the y direction, 

To maintain the same 3% lim.it, 

■OG; i 0,03 G} 

or 

± 0.09 <^y (8) 

Also in the Z direction, 

Q < 0.0 3 0} ,9) 
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The most demanding of these is the first. In order to 
model the behavior at the notch tip as a uniaxial specimen, 
only regions may be considered where the radial stress is one 
per cent, or less, of the hoop stress. 

For a circular hole in an infinite sheet, the location 
v/here G'y. = 0.0^ is at r 1. OlOl . For a two 

inch diameter hole, the distance from the edge of the hole 
where the error exceeds the stated allowable is 0.01 in. 

This precludes the use of a strain gage to make local measure 
ments. If specimen sizes were enlarged ten times until the 
uniaxial tensile region were large enough to accommodate a 
1/8" strain gage, the hole would be twenty five inches in 
diameter and the specimen would be four to six feet wide and 
ten to twelve feet long. To circumvent this difficulty, 
photoelasticity was used as the experimental method for 
measuring local strain conditions at the notch. There are 
similar size problems, but they are not quite as severe as 
for strain gages, plus the bonus of obtaining the entire 
strain field at the notch tip made photoelasticity attractive 

PHOTOELASTIC MATERIAL PROPERTIES 



SELECTION 

7075-T6 aluminum was selected for the study of notch 
fatigue since it is a popular alloy of the aircraft industry 
and one for which there is a sizeable fatigue data base. 
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To coat 7075-T6 with a birefringent material that may be 
loaded into the plastic zone, Photolastic' s PS-1 material was 
selected. It has a maximtim strain capability of 10% and is a 
good all around photoelastic material. 

Mil Handbook 5 gives typical data for 7075-T6 sheets 
of 0.080 in thickness. To extend into the yield range, strains 
of 10,000 ps must be measured with excursions up to 20,000 ps 
probably anticipated. From equation (A-1) of Appendix A, it 
is readily seen that 



If something like five fringes are desired at the 10,000 ps 
level, the value needed is 0.002. This would allow 10 

fringes at the maximum strain level expected. 

From equation (A-4) , the material fringe value is given 

by 



For PS-1, the company supplied a value of K = 0.15, and for 




£ 

M 



( 10 ) 




( 11 ) 



white light the wave length is = 22.7 x 10 ^ in; thus 

-4 

the material fringe value is 1.5 x 10 . Using equation (A-1) , 

the thickness of the coating desired is 



W<. - il - ^ = o.osa 

2 Ft 
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The coating material PS-lC was purchased in sheets nominally 
.040 in thick. The bonding agent used was that suggested by 
Photolastic. At first PC-1 was used, but for higher 
elongations it was necessary to go to PC-8, even though the 
first cement was supposed to have the elongation capability 
required. 

CHARACTERIZATION OF PS-lC PHOTOELASTIC MATEPJAL 

A careful characterization of PS-lC was conducted by 
Stenstrom (1) using two uniaxial tensile test specimens 
prepared from a single sheet of the photoelastic material. 

The first specimen was 1.0 inch wide and was used for thick- 
ness measurement and longitudinal strain data. The second 
2.0 inch wide specim.en was used for transverse strain data 
and for measuring the photoelastic coating stress fringe 
value, . 

Thickness Determination 

One longitudinal edge of the 1.0 inch wide specimen was 
lightly sanded to create a smooth, uniform cross-section of 
the photoelastic material. It was then mounted such that a 
traveling microscope could be used to measure thickness. With 
the contrast in appearance between the reflective coating and 
the base material, it was possible to accurately determine 
the thickness of each. For the specimen used, the base 
thickness was measured to be 0.0382 inches and the reflective 
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coating to be 0.0051 inches. These numbers were obtained by 
measuring eight separate points along the edge and averaging 
the results. The total average thickness of 0.0433 inches 
fell within the range of the manufacturer's measurement of 
0.044±0.002 inches. It was therefore concluded that the 
labeled thickness included both the base and the reflective 
coating after two similar specimens of different marked 
thicknesses were checked. For calculation purposes, the 
thickness of the coating was taJcen to be 0.0 382 inches since 
the path of polarized light does not penetrate the reflective 
coating of the photoelastic material. Inherent in using 
this thickness is the assum.ption that the reflective coating 
does not contribute to the reinforcement of the base material. 
Young's Modulus 

Two lines were lightly etched on the surface of the 1.0 
inch wide specimen in the horizontal direction in order to 
accommodate a Tuckerman Optical Strain Gage. The specimen 
was then loaded in tension. The readings were taken with a 
autocollimater and recorded. Using the American Instrument 
Company calibration factors, the strain was calculated with 
the following formula taken from Ref. 2. 



P • A- 

\ooo e 



( 12 ) 
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where S = Strain in microinches per inch 

F = Lozenge calibration factor, = 1.002 
A = Auto-collimater calibration factor, = 1.004 
L = Lozenge length, 2.0 inches 
R = Net reading 
e = Gage length, 1.0 inch 

With the measured strain, the measured load at each point 
and the known cross-sectional area, Young's Modulus, E, was 
determined from the slope of the stress-strain curve shown in 
Figure 1, using a linear regression analysis. 

(Tr SST.oT'i.iii <“) 

6,1 0.'598S 

where R = correlation coefficient. 

Dropping the non-zero intercept. Young's modulus was taken 
to be 

521,079 psl ( 14 ) 

Photolastic Inc. information lists this modulus to be nomin- 
ally 360,000 psi . Using thin specimens for this measure- 
ment requires very accurate thickness measurements because 
E is inversely proportional to thickness. Small errors in 
thickness can turn out to produce large percentage errors in 
area and modulus. 
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PS-IC CHARACTERIZATION 
STRESS-STRAIN CUR'/E 




Figure 1 
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Poisson's Ratio 



Using the 2.0 inch wide PS-lC specinen, two lines v/ere 
lightly etched into the surface in the longitudinal direction 
one inch apart. With the Tuckerman Optical Strain Gage 
mounted in the transverse direction as shown in Figure 2 , 
the specimen was loaded in tension. The data points from 
the test were used in a linear regression analysis with the 
following results; 

(J- Z . OZ, - l> 0O1O& (15) 

0.99991 

Poisson's ratio was then determined by dropping the non-zero 
intercept and dividing the negative of the slope by E as 
indicated in Figure 1. The measured Poisson’s ratio is 

r O. S795 

which compares with a nominal value of 0.38 supplied by 
Photolastic Inc. These tests were performed in accordance 
with ASTM standards for uniaxial properties (3) . 

Coating Fringe Value in Terms of Stress, F^ 

Prior to etching lines on the 2.0 inch v/ide specimen, it 
was loaded in tension in the elastic region, and the fringe 
number, N, was determined using a reflection polariscope. 

With several load points and known dimensions, the stress 
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was calculated and plotted against the fringe number. Using 
the data from Table I, a linear regression was made with the 
data which yielded 

CT- •+ 1.9114 ( 16 ) 

.9993 

Dropping the non-zero intercept, the stress fringe value was 
taken to be 

F- P”^**/fy'u^C)e, (17) 

Coating Fringe Value in Terms of Strain, F^. 

With the tensile modulus and the Poisson's Ratio of 
the coating material already determined, the strain coating 
fringe value can be calculated readily. From Appendix A, 

Ft “ ^ F<r = i"IG4 5^1 o"** 

CHARACTERIZATION OF 70 75-T6 
The elastic and plastic material properties of the 
aluminum panels were established by KAISER (4) using tensile 
tests of uniaxial specimens made from the same mill run as 
the sheets to be tested. The specimens were manufactured 
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TABLE I 



PS- 1C DETERMINATION 



u 


Compensator 




psi 


Reading 


N 


0.0 


0.0 


0.0 


202.37 


23.0 


0.489 


408.42 


39.0 


0.830 


610.26 


63.0 


1.340 


785.26 


79.5 


1.691 


990.53 


102.0 


2.170 


1189.21 


123.0 


2.617 



LINEAR REGRESSION; 

0 - 456.6661N + 1.92740 
R = .99932 
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and tested according to ASTM standards E8 and B5 57. ( 3) The 

specimens were cut transverse to the rolling direction, and 
the tests were performed in a 100 Kip MTS machine. The 
duration of the tests was three minutes. MICRO-MEASUREMENTS, 

E A-13-125AD-120 , precision strain gages with a temperature 
compensated bridge circuit were used on all specimens. Trans- 
verse gage sensitivity errors were corrected according to the 
manufacturer's recommendations (5). Critical cross-section 
measurements were made with a micrometer. 

YOUNG'S MODULUS (E) 

Tests were conducted using the specimen shown in 
Fig. 3 on the MTS test machine with 10,000 lb INSTRON grips, 
which gripped the specimen evenly. The results of testing 
three specimens are shown in Tables II to IV where all strain 
values have been corrected for transverse sensitivity. 

Linear regression of the first eight points in the elastic 
range of all three specimens gave 

d = ^ '-^“7 X 10^ t. (18) 

2 

with a correlation coefficient of R = 0.9997. Ignoring 
the non-zero intercept as specified by ASTM 

(19) 
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FIGURE 3 Characterization Specimen 
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TABLE II 



MTS SPECIMEN A TEST RESULTS 



2 

Cross-section = 0,03975 in 
Load, lbs, Strain, e,10 ^ in/in 



256 


615 


503 


120 4 


750 


1801 


1005 


2423 


1255 


3042 


1505 


3665 


1778 


4352 


2003 


4925 


2252 


5587 


2508 


6355 


2755 


7365 


2905 


82 30 


2984 


9045 


3037 


10150 
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TABLE III 



MTS SPECIMEN B TEST PESULTS 



2 

Cross-section = 0.03975 in 
Load, lbs. Strain , e , 10 



500 


1250 


750 


1860 


1000 


2450 


1250 


3080 


1500 


3700 


1750 


4320 


2000 


4950 


2250 


5630 


2500 


6380 


2750 


7340 


2900 


8200 


3000 


9250 


3050 


10200 


3100 


11550 


3125 


12 800 



in/in 
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TABLE IV. MTS SPECIMEN C TEST RESULTS 



Cross - 


section = 


0.03975 in 


Load, 


lbs . 


Strain 


2 72 




663 


503 




1220 


76-2 




1850 


100 8 




2 451 


12 31 




3005 


1506 




3690 


1755 




4318 


2000 




49 40 


2255 




5607 


2503 




6325 


2750 




72 30 


2900 




8120 


3000 




9200 


3060 




10500 


3100 




11500 


3125 




12500 



, 10 ^ in/in 
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POISSON'S RATIO (v) 



Tests were conducted using the specimen shov/n in 
Fig. 4 on the RIEHLE test machine with 10,000 lb RIEHLE 
grips. The results were tabulated in Table V. The trans- 
verse strain measurements were corrected by the manufacturer's 
data on transverse sensitivity. Both strains were read at 
the same load and plotted versus one another. A linear re- 
gression was made on the data resulting in; 

^ r 0.0000 — O.^ZSU I (20) 

with a correlation coefficient of R = 0.99998. Disregarding 
the constant factor as insignificant to four places 

r 0. ilSfc (2i) 

YIELD STRESS AND STRAIN HARDENING MODULUS 

The values of yield stress and strain hardening modulus 
v/ere determined graphically using the data from Tables III 
and IV. Plastic region data in Table II is not reliable 
because of excessive creep encountered during that test. 

0.2% offset yield stress, = 76,000 psi 
strain hardening modulus, E^ = 566,000 psi 
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FIGURE 4 



2 GAGE SPECIMEN 




I 
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LOAD 

lbs. 

500 

1000 

1500 

2000 

2500 

3000 

3500 

4000 

4500 

5000 

5500 

6000 

6500 

7000 

7500 

8000 

8500 

9000 

9500 



V. REIHLE SPECIMEN TEST RESULTS 



2 

Cross-section = 0.12 in 



STRAIN, £3_ 
10 ^ in/in 


STRAIN, £2 

10 ^in/in 


345 


-115 


720 


-240 


1114 


-364 


1501 


-492 


1895 


-614 


2291 


-745 


2705 


-871 


3095 


-998 


3515 


-1126 


3912 


-1255 


4332 


-1384 


4750 


-1520 


5188 


-1647 


5595 


-1784 


6075 


-1928 


6649 


-2108 


7385 


-2375 


8663 


-2888 


12245 


-4435 



23 



The graphical fit of these values to the test data can be 
seen in Fig. 5. 

RAMBERG-OSGOOD COEFFICIENTS 

The Ramberg-Osgood equation for elastic-plastic 
stress-strain characterization is given by; 




where : 

S. = strain 
(T = stress 
E = Young's modulus. 

The S and n coefficients were determined graphically from 
the data of Table IV, by the method given by Ramberg and 
Osgood ( 6) . 



■■ £ 



T L . 



o-v 



(23) 



Y\ = 



ioj (’Vt) 

\ <r 



(24) 



Lines corresponding to slopes of 0 . 7E and 0.85E were 
laid out on the X-Y recorder plots of load versus longitudinal 
strain. Load values where the rays intersected the load- 



24 



5TRC55, 



FI EURE 5 

707E-TS ALUMINUM 5TRE55-5TRA ! N CURVE 



32 j 1 1 1 1 

' E.“566»000psi — y 

\ T 
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strain curve were 



p = 3060 lbs P = 2 950 ^bs 



2 

The cross-sectional area from Fig. 3 is 0.03975 in . This 
gave 



7 7,484 psi. 




14 , 2.14 



The Ramberg-Osgood Parameters are given by 

^ = 1.479 10^^ 

n-- 

The stress-strain law for uniaxial tension throughout the 
entire range is 

-8 ^ . ,.rto3 Xl.SS 

9.86IK10 C <J- 



This is shown plotted in Fig. 5 along with data points from 
two of the specimens. 



REINFORCEMENT 

When specimens are made from thin sheets of material 
to duplicate wing skins, the coating material can provide 
some reinforcement to the specimen and influence the stress 
values. This is usually small, but since it is calculable 
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by classical developments, allowance can be readily made 
for it. 



Consider an infinitesimal elem.ent in a biaxial stress 
field oriented along principal stress directions as illus- 
trated in Fig. 6 . If the total loads applied to the 
element are and P2, they can be equated to the loads in 

the coating and the specimen. 



If there were no coating, the stress produced by P^ and P2 
would be 



The stresses in the ccm.ponent parts of the coated specim.en 
are given by 



Equating the loads in the coated and the uncoated elements. 





( 26 ) 





( 27 ) 





( 30 ) 



( 31 ) 
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FIGUEE 6. Reinforcement Model 



28 



Simplifying terms, 



U. " 






(T.c 


(32) 




Gti + 


h± 

I 


(j\c- 


(33) 













Subtracting the two equations 

cT. -- 



iTA- 



(34) 






From the plane stress constitutive equations. 






(35) 



equation (34) can be transformed into expressions for strain 



1^0^ 14'^S Kj \+0>c. 

Recalling the fundamental characteristic of the photoelastic 
coating method. 

(37) 



The reinforcement equation simplifies to 



t I 4- (ixc." 



(38) 



29 



The quantities in the bracket have all been measured; there- 
fore all photoelastic measurements must be multiplied by 
this factor to obtain uncoated specimen response 




( 39 ) 



POISSON’S RATIO MISMATCH 



As mentioned previously the notch tip may be considered, 
under certain restrictions, as a uniaxial specimen. As 
the tip is strained, the lateral contraction is different in 
the specimen and in the coating. At the interface of the 
coating and specimen. 



Thus there is a transverse strain gradient through the thick- 
ness of the coating. The birefringence will exhibit the 
average strain readings, which will lie somewhere between the 
values at the interface and the values at the surface. At 
the interface. 






( 40 ) 



at the free surface of the coating 



tec. - 



( 41 ) 
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tu. - ^.c - F^N 



(42) 



At the free surface. 






- €.tc ' 



M 



(43) 



The actual fringe value, N, will be an average bounded by 
these two, 

TT ^1. 



(44) 



Experiments with uniaxial tensile specimens have shown (7) 
that the strain variation through the thickness are only 
present near free boundaries in the plane of the specimen, 
such as the notch surface. This arises because of the edge 
effect . 




Describing the displacement field sketched above, we assume 
it to be of the following form: 
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U-- K i) 



M = V‘(y) 



UJ = UJ (i) (45) 



W 

At X = /2 there exists a contour for u which describes a 

curve in z 






2 2 

For instance g(z) = (z - h^ ) yields a good solution; 
however for now we take a more general approach and only 
assume that g(-h ) =0. 

This contour is continuously modified as one proceeds to the 
inner regions. This will be represented by 

( 47 ) 



Since the displacement in the coating is not zero at the 
specimen boundary but rather produces a constant strain equal 
to that of the specimen, another term must be added. 

utx.i) : 4W^Ci) + (48) 



The other displacements are chosen to give constant strains 



v(y) - C^y 


(49) 


lO ( ^ 


(50) 
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The strains become 



£,^z ^ ^ -TwaU) + c, 



- 

^yi 



}j£ 

^vc 

TT 



- C: 



}w . Cj 



- , }\T _ o 

' Jf 7^' 

r 4 ^tu - -f (x) a '( i) 

- 4 - O 

71 77 



(51) 



The stresses become 



£c 

^ ■ itO^ 



o:-- 



f'g +Ci i- \ 'g + + 



G*\j - 



£c 



l+Oc 



Cr 


+ 


-f + C, + Ct. + ^5 


mm ^ 


i i 





(52) 



(T, - 



£c 



^ i + > 3 < 



. \ 






'<3 + C, + 

J 
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^xe. 6 + g 



T-yi -- O 



7 - - O 

%'j 



The general form of the stress expressions is; 



(jr. ~ PW6(i) + C°>'iUr\i 

'•i 



(53) 



Using this general form as a guide, a more general solution 
is sought using the field equations. Assume txy = ^xz “ 
and y =0. The Equilibrium Equations simplify to 









7T 


4 




^>X. 


Cli 



(54) 



These equations can be solved with a stress function by 
prescribing that 

. - cr. r 



(T- 

Jv- 






^ X c> i 



i - (55) 
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rhe Compatibility Equations reduce to 

^ e , o 



vX + 


1 




O 




1 

uOc. 




1 

1 40c 






= O 



'3^9 



= O 



e= g; vcr. +cTj 



(56) 



Adding the first three compatibility equations gives 



( 



I + Oc. 

Since 9’^CTy =0 






(57) 



Then 



+ cr^) = o 



(58) 



Substituting the stress functions into (58) , we obtain the 
governing equation for it. 

|> = O (< 

Using the biharmonic theorem, 
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( 60 ) 



where 



vy = 



and 






Using separation of variables for and ^ 



f{x)g(£j 



/= K(x)k(£) 



(61) 



the stress function can be represented as 



§(x,£)= £f(X)g(i) 4- K(^)k(i) 



(62) 



The stress of principal concern here is . 



a; - -f(x) dl[£3(£)j + kw k"(z.J 



(63) 



To make it match the general form of equation (53) , let 



f(»)= I 



. c, 



cl t 



(64) 



g(z) is found by simple integration and h(x) and k(z) are 
found from solving LaPlace's Equation: 

lr\(x^ - A, Cos>K Ax + A2 Stl^k Ax 



lc(2.) = ^1 Cos -^2 : +62 6>CK\. a £ 



(65) 

(66) 
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where is the separation parameter. The stress function 

is now determined. 

^ +. Cj 2; + C 3 + (a,Co^ Ax + Aj^St-rsW 

2 . 

Differentiating, the stresses becom.e 

(J’^ : C, - Cobh Xx + Sink Cos A-2 +• S Ia. 

5 '.-Ak >/)(£> t Cob Ai -b 5 ^SLa_ ^i) (68) 

Sb^kAx + Az_ CoskAx)(75, S‘u\ -r 5^^ Cos At) 



The normal stresses must be symm.etric in therefore, 

A2 = 0 

Without loss of generality, Aj_ = 1- This simplifies the 
stress expressions somewhat to give, 

(jl - C ■" C c?bk.A^ 'A (^B| Cos A £ + £>1, S (.A. As) 

(j^ : Cosk kx (6>, Cob Ai -t- Si. A. As) 

T. , - A'" Slnk^x (b, SLa. - 6>2,CobA^:) 

>*-c 
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We require that the shear stress be zero on the free surface 
z = 0. To accommodate this,B 2 is set equal to zero. The 
stresses simplify further to the form 

(T^ - C., - Ccs XiE 

(j^ - fe, C,os.k Ax. Ai ( 70 ) 

T- - A’- S>, it-AAx 

I ' '-’I 

Going back to the compatibility equations, others can now be 
satisfied. To start with, 

O - Cu 4 <Tj =. Cl + (Ty (71) 

Because V =r O from equation (70), equation (56) and 
(71) give 

^ =0 (72) 

2>-C 

Integrating, 

(fy r X PCi) + 6(^r) (73) 
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Substituting into equation (56) 



X f"{i) -w Q"(i) = o 



( 74 ) 



The only way that this can be satisfied for arbitrary values 
of X and z is for 



F'‘(7r) = 0 and 6“(2.)=<=> (75) 



Then upon integration, 



Cl(i) - C-, ,76) 



and 



(J^j - X (Ca 2: ■+'Cs) 4- 2. t Cl) 



(77) 



Since - O from equation (70); then equation (56) 

gives 

te C fc, + (C^H-Cj) -o (78) 



This is satisfied for all values of the constants. Since 



- O from equation (56) 

^ ^ - "2 (2., 4- Cs) -^'((^^^'‘'^ 7 ) 



rO (79) 
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This equation is satisfied if = 0; therefore, all of the 
field equations have been satisfied and 



(jy- Cs-X + + C] 



( 80 ) 



The strains in the plane of the coating are now of the form 









C,-6,(i+’^c)(2o5K>^x (c.sx fCfc-z,+-C7) 



(81) 






£ 4- C 7 



The loading boundary conditions are; 2. - 



This gives 



(S 



To satisfy these requirements, we set 



At JL 



CLi\£i C^-O 



(83) 
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Cg is set to zero to seek a solution for a constant cy. 
Since =-'0s / equations (82) simplify to 



- 'Of, ^y5 - 






C I " "l)c C 1 

C^ 



( 84 ) 



Solving for the constants 



(85) 

Cl ^ (i- 



The last constant is found by setting the average trans- 
verse stress equal to zero at the boundary of the specimen. 
Assuming the specimen to be W units v;ide , (5^ = • 



From which. 



(86) 



a, - E-c^ys 5. CcA]^ -o 



Ikt 






zkl £o ^ys 



IK Coi>K„ 



( 87 ) 
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The stress expressions in final form are; 







i-E 

1 



C-ei 2^ 

c;xt|g) 




(t- OcA) «-ys 

6- V) 



( 88 ) 



a and x are neglected as being approximately zero 
compared to g^. As can be seen from the expression, if 
V = V , the stress state becomes uniaxial in the limiting 
condition; however, if there is a Poisson's Ratio mismatch, 
and there normally is, a transverse stress is created in the 
coating which influences to some degree the accuracy of 
the readings. This can be accounted for through this analy- 
sis by formulating an expression for an equivalent Poisson's 
Ratio. To do this we first formulate expressions for the 
average strains through the thickness of the coating, which 
will be what the fringes will be measuring. 







I) 




["p3E) 














( 89 ) 
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An equivalent Poisson’s Ratio is defined in the obvious way. 






i) ^ ^ 



h 



(90) 



Substituting elastic material properties referred to earlier 



i)-- 



.OS^l 



\.Oo^4r + .OZ.04 



Cot>lv (^c) 



C ® t>k C 2.Uc) 



(91) 



A plot of this expression was made assuming that W>>h^; for 
instance, let W = 20h^. A plot of the behavior is shown in 
Fig, 7. 

If the specimen material is taken into the plastic range where 
- 0.500, the expression for Poisson's Ratio becomes 



47^1 - . 1 101 



l) ^ 



Co-bk. L z.kci 



(92) 



. oA^2> 



C® 



For the same parameters of VI = 20h^, the plot of equivalent 
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Poisson's Ratio for the plastic case is given in Fig. 8. 

In both cases the combined specimen has at the free boundary 
the Poisson's Ratio of the coating; however, this falls off 
rapidly as one moves to the interior. Within three or four 
coating thicknesses of the edge the Poisson's Ratio has 
approached a constant value, but it is different than the 
specimen value alone. 

EXPERIMENTAL COMPAPJSON OF STPLSS CONCENTRATION 
IN A STRIP WEAKENED BY A CIRCULAR HOLE 
To test the validity and accuracy of utilizing the 
photoelastic method, a test case using a circular hole speci- 
men of finite width was selected. R. C. J. Howland (8) 
determined the stress distribution around the circular hole 
and with the known theoretical solution, it was possible to 
determine the accuracy of the photoelastic method. 

EXPERIMENTAL PROCEDURE 

A specimen was designed to accommodate a 2.0 inch diameter 
hole in the center. To determine a reasonable width and 
length which would allow a uniaxial field to exist between 
the center hole and the end grips, a 1/8 scale PS-lC model 
was tested under various tensile load conditions. It was 
determined that a minimum of a five-to-one ratio for length- 
to-width must exist to eliminate interference of fringe 
patterns. Fig. 9 shows the geometry selected. 



45 



.4800 



O 



•H C/3 
U 1) 

w 

cn 

• 03 




u 

jC 



0 

£ 

fO 




0 



o 



o 

o 

(0 



o 

o 



o 


O 


o 


o 


o 


o 


CvJ 


o 


00 


<• 


’il' 


lO 






^f6 



CIRCULAR HOLE SPECIMEN 



T 



40 

inches 







Data 



t = 



Points 



0.0916 in. 



Figure 9 



47 



The aluminum specimen was placed under an 8000 pound 
load in the ffCS machine, and com.pensator readings were taken 
with a Model 032 Photolastic Compensator at the 90° position 
on each side of the hole. They were found to be 122 / 121 
counts, or an average of 121.5 counts. The strain-optic 
equation from Appendix A can be written as 



£y - ti. ~ Rj- 0 M 

~eT 

At the circular boundary, - - l)c 



(93) 



X - f'<r ^ 

ty - _ 






(94) 



Making the correction for reinforcement, 

- o. 00554 « 



(95) 



The local hoop stress is 



£5 £ 



tv. 



j>9c 



(96) 



Howland's results for this case gave 34,719 psi at the edge 
of the hole. This indicates a difference of 2.4%. Later 
we will show that this difference is reduced by more than 
a factor of three by extropolating to determine the surface 
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stress from readings in the interior. 

OBLIQUE INCIDENCE 

In order to determine stress values away from the edge 
of the hole, oblique incidence was used. Fig. 10 is a 
schematic showing the coordinate system for rotation about 
the Y axis. The governing equations become, 






(ty 




(97) 



(98) 



y 



i 



Figure 10. Oblique Incidence Coordinates 
for Rotation about the y axis. 
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The strain in the x' direction is given from the strain 
transformation laws 
/ 



6.)C — Co^ 0 + til © 



( 99 ) 



For plane stress, the following stress-strain relationships 
hold: 



(die 



( 100 ) 






£ c 



l-Oc 



t, + U 



( 101 ) 



CTyc - 



£c 

7^" 






( 102 ) 



Putting these together, 



6 ^- - 
^ l-Oc 



(103) 



and the expression of strain along the oblique axis becomes 



tl-- 

V -Oc. 



(104) 
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From Fig. 10 it can be seen that r- K. . Putting 

all of these expressions into Equation (98) , an equation for 
N in terms of known quantities can be obtained. 



M 




■ 



{\ - i)c Cos'-e) + t>c (^c ' Co^''<9)j 



( 105 ) 



Solving this equation simultaneously with equation (97), 



- 



Ft 









n 1 - i)c Co s^©)/ 



(106) 



f z - 

^ (l +Oc) Si^v© 



NJ 






(i-o^)Ccse - 



(107) 



For rotation about the other axis, 




Fig. 11 shows a schematic of this oblique incidence 
arrangement. The strain in the oblique direction turns out 
to be of the same form as before. 



ty - ty Co^ & ^ 



tz S Q 



(109) 
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Figure 11. Oblique Incidence Coordinates for 
Rotation about the x axis. 



Using the expresion for 64 previously derived and the 
expression for K / the expression for the oblique fringe 
value becomes 



N 



- ihi 

■ ft 



Cos6_ 



£y (Coi'e - l)<) ' (1 - Cos's) 



( 110 ) 



Solving this simultaneously with the normal incidence 
equation (97), the strain expressions are: 



£ = 



K/ (Cos''©''l>c) - Mqx (l - Oc) Cos 9 



( 111 ) 
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( 112 ) 




M(i-OcC»i'©) - r49,,(i-'0.)Coi© 






Calibrating the Polariscope 

In oblique incidence the determination of the angle 
can be tedious because of changes in it at the surface due 
to the index of refraction being different in the coating 
than it is in air. The light ray is bent as shown in 
Fig. 12. It is more reliable to obtain the value experimen- 
tally. For instance, equations (106) , (107) (111) and (112) 

can be written as 




(113) 




(114) 



(115) 




A 



(116) 



Subtracting (114) from (113) or (115) from (116) , 




( 117 ) 



Since C - 3 = 1 from the basic equation (A-1) , 



C = B + 1 
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Figure 12. Definition of Angle of Refraction 
for the Polariscope 
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Solving two of the other equations simultaneously for A and B 



(ty - 

£xj hJ©x. + ^)c N^y 

-J- 



( 118 ) 



( 119 ) 



Using a uniaxial test to make the calibration tests, the 
expression for A, 3 and C become 





( 1 + oj (Mey - ^i©■>c) 



( 120 ) 



\b = 



- ^) Nay 
(l + o)(Me;» - 



( 121 ) 



C - ‘ 



( 122 ) 



Several readings were made of normal incidence and both 

oblique incidence fringes on a uniaxial specimen of coating 

material at different load levels. N. and were 

0 X 0y 

plotted versus N and a linear regression was made of the data 
Nex - 0-S\o ^ O. (123) 

z 0. 9969 
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-2,155 1.154 M 

R. -- 0.9997 



(12 4) 



Substituting into equations (120) to (122) and using 
■0 = l)<^= 0.3793, the constants were found to be 



A = 1.372 



B = 0.899 



C = 1.899 



Generalizing for Any Material 

In order to use the equations of oblique incidence for 
measurements on materials of any Poisson's Ratio, the ex- 
pressions for the constants were used to solve for the overall 
incidence angle ©. 



will be used in all further calculations. 

The expressions for strain in terms of N and N. become 




(125) 



Solving this equation, 0 = 31.85°. The other two expressions 
for B and C both gave 0 = 31.76°. The good agreement gave 
confidence in the numbers, and the average value of 31.8° 
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PHOTOELASTIC DATA FOR THE HOLE 

Normal and oblique compensator readings were taken at 
0.0625 inch intervals at the 90° position on each side of 
the hole as shown in Fig. 9 for a loading of 12,000 lbs., 
or a far-field stress of 16,124 psi. Six different compen- 
sator readings were made by two observers; i.e. each observer 
made three readings, and these are listed in Table VI. The 
average and standard deviation are also listed. 



TABLE VI 



CIRCULAR HOLE COMPENSATOR READINGS 
NORMAL INCIDENCE 



X AVG. STD DEV 



0.0625 


148 


151 


149 


149 


154 


14 8 


149 .8 


2.32 


0. 1250 


125 


130 


126 


124 


130 


125 


126. 7 


2.67 


0. 1875 


10 8 


110 


10 8 


108 


108 


10 7 


10 8.2 


.98 


0.2500 


94 


96 


97 


98 


97 


97 


96.5 


1.38 









OBLIQUE 


INCIDENCE 








X 














AVG. 


STD DEV 


0.0625 


126 


125 


124 


125 


12 8 


124 


125. 3 


1.51 


0. 1250 


102 


102 


100 


10 2 


10 4 


101 


101.8 


1.34 


0. 1875 


85 


85 


84 


85 


87 


86 


85.3 


1.03 


0.2500 


74 


75 


75 


75 


76 


75 


75.0 


.63 



By Chauvenet's Criterion, data will be excluded where 
probability of occurrence would be less than ^ . For six 
readings this is .0833. Assuming a Gaussian distribution 
of errors on the readings, this represents 1.73a away from 
the mean before the dat\:im becomes suspect. Table VII shows 
the upper and lower limits of acceptable data. 
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TABLE VII 



LIMITS OF ACCEPTABLE DATA 

NORMAL INCIDENCE OBLIQUE INCIDENCE 

X UPPER BOLHD LOWER BOUTID UPPER BOUND LOWER BOUND 



0.0625 


153.8 


145.8 


127.9 


122.7 


0.1250 


131. 3 


122.0 


104.1 


99.5 


0.1875 


109^9 


106.5 


87.1 


83.5 


0.2500 


98.9 


94.1 


76.1 


73.9 


Data falling outside 


of these 


ranges were eliminated 


and 



new averages computed. These results are listed in Table VIII. 



TABLE VIII 

SMOOTHED CIRCLT.AR HOLE COMPENSATOR DATA 



DATA POINT 


COMPENSATOR 


READINGS 


(in. ) 


Cn 


Sx 


0. 06 25 


149.0 


124.8 


0.125 


126.7 


101.8 


0.1875 


107.8 


85.3 


0.250 


97.0 


75.0 



Extrapolating to the Edge of the Hole 

Because it is difficult to measure the fringe count 
right at the edge of the hole, these data were plotted 
and extrapolated back to that point. It was found that 
log log values of the compensator readings plotted linearly 
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with X. The best fit was obtained by using the first three 



points. 

jL^ r 0.5390 - C-.2<2>S5X ( 130 ) 

with R = -.9996. Using the value for x = 0, 

Sx = ^^2.31 

For this Compensator, 47 counts constituted 1 fringe. 

N„ = 3.24 
9x 

Similarly for the normal incidence fringes, 

with R = .99996. Using the value for x = 0, 

C., = 176.81 
N 

and 

N = 3 . 76 

At the edge of the hole the radial stress, a^, is zero, 
which requires for the two dimensional plane stress geometry 
that 

Substituting equations (128) and (129) into equation (132), 
a relationship is found between and N for all at 
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the hole boundary. 



Ns, = O. S4 9<» Kl 



(133) 



After extrapolation the measured values of N = 3.76 and 
N-^ = 3.24 check out to v;ithin 1.3% of the theoretical values 

y X 

from equation (133) . 

Accounting for Poisson’s Ratio Mismatch 

These fringe values are extrapolated from interior 
regions where "i? = .3256; however, at the edge of the hole, 

O = .3793; therefore, they must be adjusted to the appro- 
priate Poisson's Ratio value. This is done by knowing that 
at the boundary the hoop stress is given by 




I 



(134) 



For a given Cg existing at the edge of the hole, the 
relation between the photoelastic readings for a different 
Poisson's Ratio is given from (134) to be 



h4> _ Ni. 

I 1 + 



(135) 



This, of course, assumes E and F are held constant. Thus 

S £ 

for our case, i), = .379 3 and 0j_,= .3256, and 



N 



_ [ . 040b N 
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The actual normal incidence fringe value at the edge of 
the hole is now 



N = = 3.91 



To find the corresponding value for the oblique 
incidence fringe, equation (133) is used in conjunction 
with equation (135) . 



M 






- 1.D405 K 



©X. 









(136) 



From this the value of Np. at the edge of the hole is 

u X 

(|.O40S)(3.i4) = 



The last correction to the data is for reinforcement, 
which for these materials is a factor of 1.0162 from 
equation (39) . 



N = 3.98 = 3.43 

9x 

Check against Numerical Solution 

A check can now be made against the numerical solution, 
which for this geometry and loading is given by a finite 
element solution developed by Kaiser (4) to be a = 51,678 psi. 
Using N = 3.98, the extrapolated photoelastic measurements 
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give from equation (134). 

^ _ s-i^jiipsL 

® 1.^195 

This value is within 0.3% of the finite element value. 

This is surprisingly good agreement. 

Check of Stress Distribut ion 

A check of the accuracy of the photoelastic work was 
also made on the stress distribution away from the hole by 
means of a finite element solution. The criteria established 
to determine uniform boundary stress distribution was uni- 
formity in nodal displacements along the loaded edge as 
discussed by Segerlind (9) . In the models used for FINITE 
ELEMENT ANALYSIS (FEA) in this report, nodal displacements 
were uniform to within 0.1%, and the resulting stress distri- 
bution was uniform axially to v;ithin 0.1% at the panel ends. 

In all cases two-dimensional, eight noded, isoparametric 
elements were used. These higher order elements cannot be 
loaded in an "intuitive” manner, as discussed by Zienkiewicz 
(10). Fig. 13 shows the nodal loading required to obtain a 
uniform surface load. 

Element Meshes 

Two meshes were developed for each panel analyzed. 

A reasonable effort was made to keep element corner angles 
as close to 90° as possible to reduce the effect of element 
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distortion discussed by Hopkins and Gifford (11) . All 
meshes modeled a quarter of the actual panel by using the 
two axes of symmetry as is common practice in FEA. The 
step from course to fine element meshes was made so that 
each element in the course mesh was subdivided into four 
smaller elements of the same type. Such a mesh subdivision 
can be expected to give monotonic convergence of results, 
as shown by Cook (12) , and allow extrapolation to results 
of an infinitely fine mesh. Figures 14 and 15 illustrate 
the element meshes used in this analysis. 

Using ADINA 

Once the mesh was developed, input data was prepared 
in accordance with the ADINA user's manual (13) . This same 
set of data was used as input for PSAPl, the plot routine, to 
check for errors and provide a graphical display of the 
element mesh. After preprocessing by PSAPl, the data is 
entered into ADINA for analysis. In the case of linear 
analysis, two types of stress output may be specified, 
nodal point or Gauss integration point. Nodal point output 
can be computed for up to eight node point stresses for 
each element. Since 2x2 Gauss integration was used, four 
Gauss point stresses were computed for each element. The 
2x2 Gauss integration is recognized as the most efficient 
integration order for this type of analysis [Ref. 10, p. 284]. 
The linear analysis used an isotropic linear elastic material 
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FIGURE 14 

COURSE MESH FOR CIRCULAR HOLES 
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model (MODEL "I" in Ref. 4) which required input of tensile 
modulus and material properties. 

For static analyses ADINA uses a time function method 
to apply loads in steps. Linear analysis loading 
was accomplished in a single step to a nominal value of 
3,000 lbs load. The stress output from ADINA is a listing 
of nodal or Gauss point stresses for each element. Since 
the only area of interest in this analysis was the distri- 
bution of stresses along the reduced cross-section, no large 
post-processing program was developed or used. All final 
computations using ADINA output data were accomplished on 
a HEWLETT-PACKARD 9 830 A calculator, using short programs 
coded in BASIC. If more extensive stress distribution 
information were desired, some form of automated post- 
processing would be necessary to reduce the computational 
workload. At a minimum, nodal stress outputs by ADINA must 
be averaged to obtain unique values of stress at nodes shared 
by more than one element. 

Richardson Extrapolation 

The use of course and fine meshes allows extrapola- 
tion to an infinitely fine mesh by means of the Richardson 
extrapolation (14) where: 




( 137 ) 
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where 



= extrapolated solution 
= solution obtained with h=h 

c 

= solution obtained with h=h 

F 

= linear dimension of course element 
= linear dimension of fine element 
= 2 (for this analysis) 

The exponent m is determined by the order of the discre- 
tization error O(h^) . Since h represents the length of an 
element, the element area is represented by h . In a two 

dimensional problem such as this, O(h^) is of the order of 
2 

h , the area of an element. In the mesh refinement scheme 
used Equation (137) can be rewritten 



<Tc 

<Tp 

Kc 

iY\ 



G 







(138) 



thus 




(139) 



Equation (139) then becomes the relation to obtain extrapola- 
ted stresses from coarse and fine mesh results in a two dimien— 
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sional analysis. Better extrapolations can be obtained by 
using three or more defined meshes, but the computational 
effort increases significantly. 

Optimal Stress Locations and Local Smoothing 

It is generally accepted that the most accurate 
sampling points for stresses are the Gauss integration 
points within the element [Ref. 10, p. 281]. In this analy- 
sis, the nodal points are of the greatest interest; thus 
a technique of local smoothing must be applied to the inte- 
gration point stresses to obtain nodal stresses as reported 
by Hinton and Campbell (15) . The nodal values obtained 
must then be averaged if shared by two or more elements. 
Circular Holes in Linear Material 

The FEA results for a circular hole in a finite width 
strip were used to validate the photoelastic procedures dis- 
cussed earlier. The results of Howland (8) were compared to 
both the Gauss point smoothed results and the nodal output 
results in Fig. 16. The stress concentration factor, cr/a^, 
is referenced to the far-field stress. The smoothed results 
give the best match to the results of Howland at the edge 
of the hole, and the only significant variation between the 
two FEA methods occurs within the first 0.25 inches from the 
edge. In order to obtain the 0.25 inch stress value for the 
coarse mesh, in the Gauss point smoothed result, a midside 
node value had to be obtained by the averaging method. 



70 



FIGURE 16 

CIRCULAR HOLE \=^.2 LINEAR RESULTS 




DISTANCE FROM HDLE / inches 
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TABLE IX, 



X = 0.25 FEA RESULTS - NODAL OUTPUT 



DISTANCE 
FROM HOLE, in, 



Ooo = 4500 psi 



'0 ' 



psi 



Cf^ ,psi 



0.00 


14745.5 


329.8 


0.25 


9447.0 


2181. 1 


0.50 


7232 .6 


1599.3 


0.75 


6304.7 


1473.5 


1.00 


5771.7 


1139.6 


1.25 


5471.2 


891.6 


1.50 


5266.9 


671.0 


1.75 


5091.7 


489.1 


2.00 


4940.2 


354.5 


TABLE X. 


X = 0.25 FEA 


RESULTS - GAUSS OUTPI 


DISTANCE 


Oco = 


4500 psi 


FROM HOLE, in. 


Oq ,psi 


0^ ,psi 


0.00 


14422.8 


121.2 


0.25 


8721.2 


1927. 1 


0.50 


7317.1 


16 39.7 


0.75 


6218.6 


1452.4 


1.00 


5779.6 


1151.5 


1.25 


5453.3 


891.0 


1.50 


5267.1 


677.3 


1. 75 


5085.9 


486.9 


2.00 


4937.9 


343.8 
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The linear distribution of smoothed stresses along the 
sides of the element, (15) , appears to produce a less accurate 
result in this area of extreme stress gradient, when compared 
to ADINA’s nodal output result. This tendency was noted 
in all cases; however, the peak stress values from smoothed 
results consistently gave better correlation with other inves- 
tigators ( 16) . 

A circular hole with X=0.25 was also analyzed and 
compared to the photoelastic experimental results. The 
stress distribution is given in Table IX and Table X. The 
data for the first two locations has been normalized and 
entered into Table XI to show a comparison with photoelastic 
results similarly normalized. 



TABLE XI 

COMPARISON OF STRESS DISTRIBUTION 



Finite Element 



Photoelastic 



X 


Oy-/0„ 


0 x/ 


0 -^/ 0 00 


0 9 / 0 ’ CO 


0.00 


-.0269 


3.2051 


-.1167 


3.0780 


0.25 


.4282 


2.0993 


.4111 


2.1627 



Comparisons of the photoelastic data with the finite 
element show agreement within 3% and 4% on the average. This 
is not as good as the correlation made at the edge of the 
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0 . 



hole where we required 0 ^ 



If 0 is allowed to be 
r 



whatever the readings show; then the correlation is in the 
3% or 4% range; however as can be seen at x = 0, the finite 
element is not exact either, especially at these small x 
values . 

Sensitivity Studies 

At the edge of the hole where ”0 = .379 3, the expres- 
sions for strain from equations (128) and (129) become 



Substituting these expressions into the plane stress consti- 
tutive equations 



ty. - I.S198 jsie-K 

^ ^ K I ~ ^ .4 ^ ^ ^^ex. 






(140) 



(141) 



- Z. 08^1 Xio^/sl - 2.45^7 
(y z. 2- ^ 9 3c> \ 



(142) 



(143) 



The sensitivity of the stress is given by 




( 144 ) 
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Applying this to equations (142) and (14 3) 



<d(T,, - 2-0597 10 *- cInI - 2.4S<37 X 1 ^ci. ( 145) 

A(T = Z-3950 ^ (146) 

Selecting the maximum deviation of N and N as 1.73a, where 

9 X 

a is the standard deviation of the smoothed data 

dM = t.2z) = ± l-l\ 

dNle;.= (|.75)(t .34) 3 ii.45 



Then da can have two values depending upon how the plus and 
minus signs couple up. The worst of these are: 

dff> - i ii»97 
d(i; : ± P=-'^ 

These figures represent the maximum errors that would be 
expected due to misreading the fringe pattern, and they are 
approximately 4% on and explains how non-zero results 
for a^ arise. 
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CONCLUSIONS AND RESULTS 



From very careful experimental procedures where no 
hcindbook values for material properties were used but all 
were carefully determined, photoelastic raeasxirements checked 
theoretical results within 4%. This error included all errors 
in material properties, in operator readings and in equip- 
ment calibration. Without very careful and detailed measure- 
ments on every parameter, photoelastic results could easily 
be in error by 15% to 20%. 

With this level of accuracy established by controlled 
laboratory procedures and calibrations of equipment in place, 
we will now proceed to study the elastic and plastic behavior 
of notches where fatigue failure originate. 
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APPENDIX A 



Equations of Photoelasticity Coatings 




(A-1) 



(A-2) 



(A-3) 



(A-4) 



(A-5) 
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